Ratchet effects in two-dimensional systems with a lateral periodic potential 
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Radiation-induced ratchet electric currents have been studied theoretically in graphene with a pe- 
riodic noncentrosymmetric lateral potential. The ratchet current generated under normal incidence 
is shown to consist of two contributions, one of them being polarization-independent and propor- 
tional to the energy relaxation time, and another controlled solely by elastic scattering processes 
and sensitive to both the linear and circular polarization of radiation. Two realistic mechanisms 
of electron scattering in graphene are considered. For short-range defects, the ratchet current is 
helicity-dependent but independent of the direction of linear polarization. For the Coulomb im- 
purity scattering, the ratchet current is forbidden for the radiation linearly polarized in the plane 
perpendicular to the lateral-potential modulation direction. For comparison, the ratchet currents 
in a quantum well with a lateral superlattice are calculated at low temperatures with allowance for 
the dependence of the momentum relaxation time on the electron energy. 
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I. INTRODUCTION 

Noncentrosymmetric periodic systems being driven 
out of thermal equilibrium by a time-oscillating force, 
stochastic or deterministic, are able to transport parti- 
cles even if the force is zero on average. This directed 
transport, generally known as the ratchet effect, is rele- 
vant to different fields of natural sciences. Various kinds 
of symmetry-breaking micro- and nanometer-sized arti- 
ficial structures have been proposed and fabricated to 
model ratchets and investigate their fundamental prop- 
erties, for review see, e.g., Refs. [THU Directed motion of 
Brownian particles in water have been induced by mod- 
ulating in time a spatially periodic but asymmetric opti- 
cal potential.-^ Electronic ratchets where rectification of 
thermal fluctuations is achieved in systems with inhomo- 
geneous distribution of temperature have been studied 
theoretically in Refs. 6-8. In semiconductor nanostruc- 
tures, the ratchet effect has been demonstrated in various 
systems with asymmetric scatterer arrays based on both 
A3B5 9 12 and Si/Ge^ materials as well as with asym- 
metric lateral superlatticesP^HU The spin ratchets have 
been proposed for two-dimensional systems with sym- 
metric periodic potential and driving force but with the 
Rashba spin-orbit interaction! 16 * 17 * Recently, the ratchet 
current induced by terahertz radiation has been observed 
in semiconductor quantum-well (QW) structure with a 
one-dimensional lateral periodic potential induced either 
by etching a noncentrosymmetric grating into the sample 
cap layer^or by deposition of micropatterned metal-gate 
fingers. 19 Quantum graphene ratchets formed by asym- 
metric periodic strain with the period comparable with 
the de Broglie wavelength of free carriers have been stud- 
ied in Ref. |20l 

In the present theoretical work we consider a classical 
graphene ratchet consisting of a graphene sheet fabri- 



metal stripes 




FIG. 1: Schematic representation of the studied ratchet struc- 
ture. 



cated on highly resistive substrate and covered succes- 
sively by a thin dielectric layer and a periodic grating 
of semi-transparent metal fingers, Fig. [I] A technologi- 
cal opportunity to grow such a system is demonstrated 
in Refs. 21 22 where photodetectors with multiple inter- 
digitated metal fingers fabricated on the graphene were 
reported. To achieve symmetry breaking the superim- 
posed lateral structure may form the infinite sequence 
...ACBC ACBC ... with A, B representing metal fingers 
of the thicknesses a and b (a 7^ 6), and C, C representing 
hollows of different thicknesses c and d , Fig. [T] 

The paper is organized as follows. In Sec. [TT] we for- 
mulate the basic concept of the problem. In Sec. |III| we 
take into account the radiation-induced heating of the 
free carriers and spatial modulation of the heating, and 
deduce the polarization-independent contribution to the 
ratchet current. In Sec. 



IV we develop the Boltzmann 



^Electronic address: golub@coherent.ioffe.ru 



kinetic formalism and obtain microscopic expressions for 
the polarization-dependent ratchet currents. The results 
are discussed in Sec. [V| In Sec. [VI] the summary of the 
research is outlined. 
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II. BASIC CONCEPT 



According to the theory proposed in Refs. I18|19l in 
the studied systems, the pulsating ratchets, the electric 
current generation is based on the combined action of a 
static spatially-periodic in-plane potential 

V(x) = V cos (qx + (fv) 

and a spatially modulated electric-field amplitude 

E(x) = E [l + ft cos (qx + <Pe)] 

of the normally-incident radiation. Here q = 2ir/d with d 
being the superlattice period along x. In the considered 
pulsating ratchets, the symmetry-breaking is described 
by the phase shift (fv — Pe different from an integer 
number of 7r, and the electromotive force is proportional 
to sin (cp v — <Pe)- The structured graphene presented in 
Fig. [I] can also show ratchet effects. The lateral potential 
V(x) can arise due to the strain, tensile or compressive, 
in graphene areas located beneath the fingers A and B 
whereas the in-plane modulation of the pump radiation 
appears due to near-field effects of the THz radiation 
propagating through the grating. For i / 5, C ^ C", 
the shapes and local extrema of the periodic functions 
V(x) and E(x) are naturally shifted with respect to each 
other resulting in a difference between the phases cpv and 
tpE- For example, for b = c = 2a = 2c 7 , a crude estimate 
gives a value of 0.4 for sin (cpy — <pe)- 

The symmetry imposes restrictions on the polarization 
dependence of the ratchet currents. Graphene modulated 
by an asymmetric lateral potential V(x) has the point 
group symmetry C^ v with the C2 axis parallel to the 
modulation direction x and the mirror reflection plane 
perpendicular toy. It follows then that the net dc current 
density components j x and j y are related to components 
of the polarization unit vector e and amplitude E(x) = 
\E(x)\ of the normally incident radiation by four linearly 
independent coefficients 



j x = [xo(\e x \ 2 + \e v f) + XL (\e x \ 2 - \e y \ 2 )] &{x) 



dV 

dx 



j v = [xL(e x e* y + e* x e y ) + -yP cilc ] E 2 {a 



dV 

dx 



(1) 



Here we use the notations x and y for the in-plane coor- 
dinates, the bar denotes averaging over the coordinate x, 
and P C i rc = i(e x e* — e* x e y ) is the degree of the radiation 
circular polarization. The coefficient xo describes the 
contribution to the ratchet current insensitive to the po- 
larization state, while the remaining three coefficients de- 
scribe the linear (xl, Xl) and circular (7) ratchet effects. 
We develop the kinetic theory allowing us to derive equa- 
tions for all these coefficients. As compared to Refs. ll8|T9l 
the theory will be substantially extended to take into ac- 
count specific properties of graphene, namely, (i) the lin- 
ear, Dirac-like dispersion of electron energy leads to a 
strong energy dependence of elastic relaxation times and 



(ii) the electron gas in doped graphene is degenerate even 
at room temperature. 

Taking the electric-field and lateral-potential modula- 
tion in the above simplest form we obtain for the average 
in Eq. ^ 



E 2 (x) 



dV 

dx 



qV^hEl sin((/v - <p E ) 



(2) 



It is worth noting that the factor ft sin (cpy — in 
Eq. ([2| should be strongly sensitive to the geometry of 
the structured coating of the graphene, namely, the thick- 
nesses of stripes and heights of the metal fingers. Fur- 
thermore, bias voltages applied to the grid elements A 
and B should substantially change the lateral potential 
V(x) as well as the relative amplitude ft and phase ipE 
of the field modulation. The replacement of thin semi- 
transparent metal fingers by thick gate stripes should re- 
veal plasmonic effects similarly to photoresponse of the 
THz plasmonic broadband detectors, see Refs. I15|23l and 
references therein. 

Hereafter we consider a graphene sheet with the lateral 
potential V(x). The electron energy in each valley, K or 
K' \ is given by 



E k = hv k + V(x) 



(3) 



where is the electron speed in graphene and the two- 
dimensional wave vector k is referred to the vortex of 
the hexagonal Brillouin zone. Since in the model under 
consideration the behavior of electrons in the K or K' 
valleys is identical we consider the current generation in 
one of them and then double the result. In the course of 
presenting the results we will supplement them with simi- 
lar results obtained for the QW ratchet with the electron 
parabolic dispersion E k = ft 2 /c 2 /2m + V(x). For ade- 
quate comparison of the two low-dimensional systems we 
have extended the theory of Ref. |4] for QW ratchets to 
take into account the dependence of the momentum re- 
laxation time on the electron energy and the degenerate 
statistics. 



III. SEEBECK RATCHET CURRENT 

In this and next sections we will successively consider 
two mechanisms of the ratchet current. In the Seebeck 
ratchet effect, the spatially- modulated radiation heats 
the electron gas changing its effective temperature from 
the equilibrium value T to T(x) = f + 5T(x). Here f 
is the average electron temperature and ST(x) oscillates 
in space with the period d. In turn, the correction ST(x) 
causes the inhomogeneous correction to the conductiv- 
ity 5a(x) = (da/dT)ST(x). Bearing in mind Ohm's law 
j = <jE, replacing the dc electric field E by —(l/e)dV/dx 
and a by 5<j(x) we obtain for the ratchet current 



1 , dV(x) 
J * = V\ 5a{x) ^x- 



(4) 
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Here e < is the electron charge. The nonequilibrium 
electron temperature can be found from the energy bal- 
ance equation 



T(x) 



= hu>G(x) 



(5) 



where r £ is the electron energy relaxation time, uo is the 
radiation frequency, the temperature is expressed in en- 
ergy units, and G(x) is the photon absorption rate per 
electron. From Eqs. Q, (|5| we derive the working equa- 
tion 



hwr £ <9cr dV(x) 



\e\ dT 



dx 



(6) 



which can be used for a degenerate two-dimensional gas 
in both graphene and QWs. In graphene, the Drude 
absorption rate per particle is inversely proportional to 
the Fermi energy e F , 



2„,2 



G{x) 



e v. 



T tr 2E 2 (X) 



e F 1 + (wr tr ) 2 hw 



(7) 



while for QWs Vq/e f should be replaced by the inversed 
electron effective mass 1/ra. Here r tr is the transport 
relaxation time which determines the low-temperature 
conductivity. 

The temperature dependence of conductivity for de- 
generate electron gas in graphene is well documented^ 



do_ 
dT 



Tie 
3^2 



Ts F 



while for the QW systems one has 



da 



QW 



dT 



(8) 



(9) 



Here primes denote differentiation over the electron en- 
ergy e = £k = hvok or h 2 k 2 /(2m), and r\{e) is the mo- 
mentum relaxation time of a nonequilibrium correction 
to the electron distribution function depending as cos cpk 
on the azimuthal angle cpk of the electron wave vector k. 
Note that the transport relaxation time r tr is equal to 
ti(£f). 

From Eqs. (|6|-([9| we finally obtain the Seebeck con- 
tribution to the polarization-independent current Eq. 
described by coefficient xo- For graphene it is given by 



Xo 



2ire 3 v^Tr £ 



3h 2 s F 1 + (oJT tr ) 2 
while for the QW structures one has 



v S 

A0,QW 



2ire s Tr F 



3h 2 m 1 + (W tr ) : 



(10) 



(11) 



The analysis of Eqs. (10) and (11) for different mecha- 
nisms of electron scattering is postponed to Sec. |V| 



IV. POLARIZATION-DEPENDENT RATCHET 
CURRENTS 

In the presence of normally-incident radiation, an elec- 
tron is subjected to the periodic force 



F(x) = e [E(x)e 



■ C.C. 



dV(x) 
dx 



(12) 



where x is the unit vector in the x direction. In the pre- 
vious section we could avoid the consideration in terms of 
the Boltzmann kinetic equation for the electron distribu- 
tion function and, instead, used the known expressions 
for the electron conductivity. However, the second mech- 
anism of the ratchet currents should be treated on the 
base of the Boltzmann equation 



d_ 
dt 



d F[x) d 

dx 



h dk 



fh(x) + Qh(f)=0. (13) 



Here is the velocity of an electron with the wavevec- 
tor k equal to v^k/k in graphene and hk/m in a con- 
ventional QW, and Qk is the collision integral. In what 
follows we assume that r tr ,cj _1 <C r s and neglect the en- 
ergy relaxation in Eq. (13) in which case the integral Qk 
describes only momentum relaxation processes. Thus, in 
the second mechanism the ratchet current can be inde- 
pendent of the energy relaxation time whereas in the first 
mechanism related to the carrier heating the current is 
proportional to r £ . 

In terms of the distribution function the electric cur- 
rent density in graphene is written as 



3 = ve^2v k f k (x) 



(14) 



where the factor v accounts for the spin and valley de- 
generacy, in QW structures v = 2 and in graphene v = 4. 

In order not to overload the theory with too cumber- 
some equations we impose the following properties of the 
system under consideration: the electron mean free path 
h — ^o^tr and energy diffusion length l £ = vo ^r tr r £ are 
both small compared with the superlattice period d; and 
the ac diffusion is neglected which is valid if vq <C ud. 
On the other hand, no restrictions are imposed on the 
value of co>T tr . Moreover, we assume the radiation electric 
field and the lateral potential to be weak enough, so that 



\eE \v r tY T and 



V(x) - V{x) 



Then, according to the phenomenological Eqs. ([T]) the 
function fk(x) should be calculated in the third order 
of the perturbation theory, including the second order 
in the electric-field amplitude and the first order in the 
lateral potential. Taking into account that our aim is to 
derive expression for the sum ( 14) rather than to find the 



function fk(x) explicitly we can express this sum in the 
form 



^fi E /\x)E*(x) 



k 



dk 



C.C. 



(15) 
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where a = x,y, and fj^ V ^ is the second-order iteration 
linear both in E(x) and dV(x)/dx. It can be found from 
the equation 

-■^ + v k , x l)fi E J ) ^) + QM (EV) ) 



e E{x) df (V) , dVldfW ^ 
h dk dx h dk x 

where the first-order corrections are given by 



(16) 



ft > = V(x)&(e k ), f k »> = -e Tlu E{x)v k ae k ) . (17) 

Here fo(£k) is the equilibrium Fermi-Dirac function at 
V(x) 0, and r luj n/(l - ioun). 



A. Graphene 

In contrast to systems with the parabolic energy dis- 
persion, in graphene the relaxation time t\ is energy- 
dependent even for the short-range scattering potential. 
Moreover, the derivative dv a /dk^ is ^-dependent which 
means that the current (Jl5| cannot be in general ex- 
pressed exclusively in terms of the macroscopic fluctu- 

a tionSN u (x)=4Zf i k E J ) - 

k 

To calculate the ratchet current, we solve the kinetic 
Eq. (16) and find a contribution to fj^ V \x) even in k. 
It is convenient to present it as a sum of isotropic part 

(/Iw^W) an d anisotropic part 5fj^ V \x), where the 
angular brackets denote averaging over the directions of 
k. The isotropic part describing a nonequilibrium correc- 
tion to the energy distribution of electrons has the form 



Here = — io;, where r 2 is the relaxation time 
of the second-order harmonics of the distribution func- 
tion proportional to cos2(/?fc = (k 2 — k 2 )/k 2 or sm2cp k = 
2k x ky I k . 



Substitution of the solution f£ V) = (f^j+Sf^ 



into Eq. (115) yields the ratchet current exactly in the 
form of Eq. ( l|) where, see the details in Appendix, 



3 2 

Xo = JJr (ReSi - lmS 2 ) 



7 



27rh 2 
^(ReS 2 -ImSi) 



(20a) 
(20b) 
(20c) 



2tt/i 2 

The complex coefficients Si ^ are defined by the following 
expressions 



Si 



s 2 = 



3 / r l 



T2 - ( —) - 2 £ kt) T2uTh 



(tie)' (ti u e)' 



COS 



(21) 



where one should set e — £f- 

One can see that, for the second mechanism, the 
ratchet current reveals contributions both dependent and 
independent of the polarization. 

In order to reveal the physical nature of the 
polarization-dependent ratchet currents we will consider 
in more detail one of the contributions to these currents, 
e.g., the contribution due to the first terms in the square 
brackets in Eq. (18). Substituting this term to Eq. (15), 
using the average (A2) and replacing (— f ) by 5(s — Sf), 



see Appendix, we can reduce this contribution to the form 



AEV)\ 
Jkuj 



(~fo) 



r 

E x {x) 



iev c 
~2u 
dV 

dx 



(18) 



for lu V(x 



dE x 
dx 



In fact, Eq. ( 18 ) describes the local oscillation of the elec- 



tron kinetic energy induced by a combined action of the 
ac electric field and the dc static electron potential. The 
oscillation amplitude increases with decreasing the fre- 
quency uj until the latter becomes comparable with r" 1 
whereupon the amplitude is stabilized by the energy re- 
laxation processes. 



AEV) 



describes 



The anisotropic oscillating correction uj koJ 
the dynamic alignment of electron momenta in the k- 
space. It has the form 



(EV) _ eVQT2u 



) dV 
(E x cos 2(p k + E y sin 2<p k ) — 



(19) 



fdE, 



dE, 



+foTiujV(x) — — cos2(p k + — - sin2(p k 



V dx 



dx 



e 2 v? 



HneY^SN^E^x) + c.c. (22) 



2e 

Here SN u (x) is the electron density fluctuation re- 
lated to the first term in Eq. (18) and proportional to 

E x (x)(dV/dx): 



\pii 2 d\f 
SNM) = 2^( Ti - £ )U f 5(£f)^(x)— , 

where we introduced the density of states 



9(e) 



h 2 v. 



2 • 



(23) 



(24) 



It is convenient to hold the further interpretation by us- 
ing the illustration sketched in Fig. [2j The lateral poten- 
tial is chosen in the form of Vq cos (2irx/d), as shown 
in Fig. [2^i, and the phase of the electric-field ampli- 
tudes E XjV (x) is shifted by cpv — ^>e — k/2 as follows: 
E Xy y(x)/E x ^ y — l+/isin (27rx/d) as illustrated by Fig. [2)3. 
In this case the products E(x)(dV/dx) and 



E 2 (x) 



dV 

dx 



2EE(x) 



dV 

dx 
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E(x)/E = 1 + hsin (2nx/d) 




FIG. 2: (a) Schematic representation of the spatial variation 
of the Dirac point along the x axis due to the built-in lateral 
potential V(x) (solid curve), see Eq. (J3|. Dashed- and-dotted 
horizontal line shows the equilibrium Fermi energy, (b) Spa- 
tial variation of the first derivative of the lateral potential and 
spatial modulation of the radiation electric field. 



are nonzero. The time and space variation of SN is rep- 
resented by 

5N(x, t) oc Eq ^— (lm£2 cos out + ReS2 sin out) , 



where E = E is the electric-field scalar amplitude, and 
62 is given by the first term in the second equation (21). 
For the circularly polarized light E x (t) oc cos out, E y (t)(x 
smout and, therefore, the time average of the product 
5N(x,t)E y (t) is proportional to Re^. On the other 
hand, the xx-related current in Eq. is induced by 
the linearly polarized light with E x (t), E y (t) oc cos out re- 
sulting in xl oc ImS2, in agreement to Eqs. (21). The 



other terms in these equations are obtained in the sim- 
ilar way by taking into account the explicit expressions 

for (fiT)™^fiT- 



B. Quantum-well structures 

In this subsection we extend the theory of Refs. 0J 
I18|19l to consider the degenerate statistics of the two- 
dimensional electron gas in heterostructures and take 



into account the possible difference between the relax- 
ation times n and T2 and their energy dependence. 
Solution of the kinetic Eq. (16) shows that Eqs. Jl8|) 



and ( |19| ) are also valid for the QW structures provided Vq 
and {jiuj/sk)' are replaced by the squared Fermi velocity 
v 2 ? = 2e F /m and r^/sk, respectively. Furth e rmor e, for 
the QW structures we also obtain Eqs. (20a)-(20c) with 
Vp instead of Vq and given by 



Si,qw = st[t2cjt[ u 



2e 



(25) 



>2,QW 



(ne)' {t 1uj s)' 



oue 



2oue 



We note that for a quadratic energy dispersion = 
h 2 k 2 /(2m) and energy-independent time t\ = r tr , the 
partial derivative d(r\y a ) /dkp = 5 a phT tr /m is indepen- 
dent of fc, and Eq. ([15]) for the ratchet current takes the 
fornP 



2e 2 r tr 



Re 



m 



where ^ w (x)=2EL 

k 

to the electron density. 



(EV) 



5N u (x)E*(x) 



is the second-order correction 



V. DISCUSSION 

We calculate the ratchet current excitation spectrum 
for two types of elastic scattering actual for graphene. 
For scattering by short-range defects one has 



Tl Ttr 



T 2 



2 ' 



and Eqs. (21) yield xl = for this case. In contrast, for 



scattering by Coulomb impurities, when 



TL T tr 



e F 



r 2 3ri 



Si = 0, and xl — Xo- This means that x component of 
the current is generated in this case only by x-polarized 
radiation: j x oc \e x \ 2 . 

In Fig. [3] we plot the ratchet current frequency depen- 
dence for both types of scattering. It can be seen that 
the coefficients xo,l and 7 have complex non- monotonous 
behavior. In the static limit, ou — >• 0, xo remains finite 
while the circular ratchet current is absent, 7^0. This 
is correct because helicity-dependent effects can not be 
present for static electric field. For Coulomb scattering 
7 also tends to zero but this occurs at ou ~ r" 1 <C r^ 1 



as discussed in the paragraph below Eq. (18). 



Now we turn to the QW structures. For scattering by 
short-range defects in QWs when 

ri= r 2 = T tr , 



Eqs. (25) yield 



Xo = Xl = Xl 



-0UT tT 7, 7 : 



27rh 2 mou 1 + (our tr ) 2 
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FIG. 3: Frequency dependencies of ratchet currents in graphene for scattering by short-range defects (a) and Coulomb impu- 
rities (b). 
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FIG. 4: Frequency dependencies of ratchet currents in QW structures at two types of elastic scattering: by short-range defects 
(a) and by Coulomb impurities (b). Inset shows the frequency range where xo changes its sign. 



Again, we obtain generation of j x at radiation polariza- 
tion along the x axis only, but, in contrast to graphene, 
this takes place for short-range scattering. Comparing 
this expression with the result for Boltzmann statistics 
valid at room temperature T room P we get 

liT = 0) T room 

rs^ m 

7(^room) £F 

This estimation implies that the helicity-dependent 
ratchet current shows no remarkable variation with tem- 
perature. 

For scattering by Coulomb impurities in QWs, when 

£ 

we obtain xo 7^ XL- Figure [4] shows the ratchet current 
in QW structures. One can see that, for Coulomb scat- 
tering, polarization-dependent ratchet currents are sign- 
constant, while xo nas a maximum and changes its sign 
at ujT tr ~ 2, see inset in Fig. [4] 



The Seebeck contribution to the ratchet current is ab- 
sent in graphene for both considered types of elastic 
scattering, cf. Eq. (10). Nonzero Seebeck ratchet cur- 



rent in graphene appears, e.g., at scattering by screened 
Coulomb potential. 

The Seebeck ratchet current is also absent for short- 
range scattering in QWs, cf. Eq. ( pT[ ). In contrast, for 
scattering by Coulomb impurities X^q W is nonzero. Its 
ratio to the elastic-scattering contribution can be esti- 
mated as 



s 

Xq,QW 
X0,QW 



At the low temperature T « 4 K, the energy relaxation 
time in QWs r e ~ 1 ns 25 while the transport scattering 
time T tr ~ 1 ps. Therefore the Seebeck contribution to 
the polarization-independent ratchet current dominates 
in QW structures at low temperatures for scattering by 
the smooth Coulomb potential. 
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VI. SUMMARY 

To summarize, radiation induced electric currents in 
graphene with a spatially periodic noncentrosymmetric 
lateral potential are studied theoretically. The ratchet 
current is shown to consist of polarization-independent 
contribution and the contribution sensitive to linear 
and circular polarization of radiation. Two microscopic 
mechanisms of the polarization-independent ratchet cur- 
rent are considered and compared: the Seebeck con- 
tribution generated in the course of energy relaxation 
and the current controlled by elastic scattering pro- 
cesses. We demonstrate that the ratchet current exci- 
tation spectrum strongly depends on the type of elastic 
scattering. Two realistic mechanisms of electron scat- 
tering in graphene are analyzed. For a short-range po- 
tential, there are polarization-independent and helicity- 
dependent currents, while the linear polarization leads 
to no ratchet current. For the Coulomb scattering, the 
linearly-polarized radiation generates the ratchet current 
only for the polarization vector parallel to the lateral- 
potential modulation direction. The Seebeck ratchet cur- 
rent is shown to vanish for both types of elastic scat- 
tering in graphene. For comparison, we have analyzed 
the ratchet effect in QW structures with a lateral su- 
perlattice and degenerate electron gas and demonstrated 
the polarization-dependent effects as well as the Seebeck 
ratchet current for the Coulomb scattering. These re- 
sults show that ratchet current measurements allow one 
to identify a dominant mechanism of elastic scattering in 
graphene. 
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The summation over k in (15) is performed in two 



stages, firstly, by averaging over the directions of the 
wavevector k and, secondly, by integration over the mod- 
ulus k = \k\ or, equivalently, over the energy e = £k = 



hvok. For finding the contribution of (^fj^ V ^ 
ratchet current, it suffices to use the average 



to the 



hvl (ne)' 



Since the anisotropic correction Sfj^ V ^ is a linear func- 
tion of cos 2(fk and sin 2(^&, its contribution to the ratchet 
current j is found by using the identities 





The calculation at the second stage, integration over 
the energy, is simplified by the assumption of the degen- 
erate statistics in which case the the first derivative /q 
can be replaced by the delta-function — S(e — £p). The 

terms in fj^ V ^ proportional to the second derivative /q 
are treated by using the identity 



Appendix A: Derivation of equations for Si and & 



4j2fo(e k )F(s k ) = [g(s)F(s)]' £ __ 



(A4) 



In order to calculate explicitly the current (15) for 
graphene one can apply the identity 



dfavp) 
dk a 



n 



Tl Y kgkp 

k 2 



(?) 



(Al) 



where F(e) is an arbitrary smooth function of e and g(e) 
is the electron density of states given by Eq. (24). 
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